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The odon. larval assemblage from Rio Pinolapa (RP) in the municipality of Tepal- 
catepec, Michoacan, is described. Sampling was conducted twice in each season (8 
trips in total), and additionally some physicochemical variables of the river channel 
were recorded. Strata (shores, riffles and eddies) and seasonal variation of assemblag¬ 
es are described and compared using classical diversity measures such as Shannon’s 
diversity index, Simpson’s diversity index as a dominance measure, Margalef's rich¬ 
ness index and Pielou’s evenness index. For comparing strata and seasonal diversity 
the Renyi’s diversity profiles were used. A Cluster Analysis was performed on a Bray- 
-Curtis similarity matrix to explore the faunal relationships among year seasons and 
strata. CCA was also performed to investigate the relationships between the physi¬ 
cochemical and species abundance matrixes. As results, 28 spp. (12 Zygoptera and 
16 Anisoptera) were recorded as larvae. Most abundant species were Erpetogomphus 
elaps, Brechmorhoga praecox and Phyllogomphoides luisi. The highest number of spp. 
was registered in winter and the lowest in summer. Among strata the highest abun¬ 
dance was recorded in riffles, although the shoreline had the largest number of spp. 
The most similar assemblages were those of autumn and winter. Shore habitats were 
more heterogeneous than eddies and riffles and this could explain the larger number 
of species. The Clench’s model explains better the data. Additionally, we used the 
slope of cumulative number of spp. curve for assessing completeness of the RP list. 
CCA was significant, with pH, autumn, shoreline and riffles the most important vari¬ 
ables. This means that species variation is related to physicochemical, temporal and 
strata conditions in RP. 
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INTRODUCTION 

The estimation of biodiversity has become one of the major goals for ecolo¬ 
gists up to now (MAGURRAN, 2004). Conservation of biodiversity requires the 
knowledge of its patterns and magnitude (BASELGA & NOVOA, 2008), where 
the process of making an inventory of unknown assemblages, as well as describ¬ 
ing new species, represents the first step in understanding that biodiversity, and 
also the foundation for any later research related to biodiversity. Unfortunately, 
this kind of work has received relatively little attention by many journals, propi¬ 
tiating a great abandonment by researchers (WHEELER, 2004; WHEELER et 
al., 2004; DE CARVALHO et al., 2005). 

Different methods for biodiversity assessment have been developed as an in¬ 
dicator of ecosystem conditions, conservation goals, management and environ¬ 
mental monitoring (SPELLERBERG, 1991). Usually, the number of species 
is the more used measure of diversity (MAYR, 1992; MORENO, 2000). How¬ 
ever, diversity depends not only on the number of species but also on the relative 
abundance of them. Generally, species are distributed according to hierarchical 
abundance classes, from some very abundant to some very rare. As some species 
become more abundant and others become rarer, the lower the biodiversity of 
the community will be. Then, the conservation of biodiversity is mainly a prob¬ 
lem concerning to the ecological behavior of the rare species. 

In this manner, measuring the relative abundance of species will allow us to 
identify those species that, due to their rareness in the community, are more sen¬ 
sitive to environmental changes (MORENO, 2001) and, consequently, more ex¬ 
posed to extirpation. 

On the other hand, Odonata have become among the most used aquatic insect 
groups in ecological quality assessment today, because they are relatively large, 
abundant, widely distributed, easy to collect, the larvae are rather sedentary and 
easy to rear in laboratory, and there is limited genetic variation (HELLAWELL, 
1986). Also, they are relatively easy to identify in comparison with other groups 
(STORK, 1994), and are generally well represented in aquatic samples (HAM¬ 
MOND, 1994). Moreover, they respond quickly to environmental stress (NOSS, 
1990). Fortunately, the identification of immature stages has been developed 
in the last two decades in Mexico. However, an important gap still remains on 
the knowledge of the structure and function of Mexican Odonata assemblages, 
with relatively few studies dealing with this subject (NOVELO-GUTIERREZ 
& GONZALEZ-SORIANO, 1991; GOMEZ-ANAYA et al., 2000; NOVELO- 
-GUTIERREZ et al., 2002; ALONSO-EGUIALIS, 2004). The exploration of 
patterns in time and space of Odonata larval assemblages can supply basic data 
for future research. 

The goal of this work is to describe the structure and seasonal variation of the 
Odonata larval assemblage from the Rio Pinolapa and relate it to environmental 
factors. 
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STUDY AREA 

The Rio Pinolapa (RP) is located at (19°00.524N; 103°01.456W), municipality of Tepalcatepec, 
in Michoacan State, Mexico (Fig. 1). The sampling site is at 616 m asl. Average gradient was 0.02 
(1°08’44.75”, n = 7), minimum = 0.006 (0°20’37.57”), maximum = 0,042 (2°24’ 18.03”). Average depth 
= 0.11 m (IC = 0.05 - 0.19 m, n = 8); average width = 2.18 m (1C = 0.65 - 3.71 m, n = 8); current ve¬ 
locity = 37.58 m/s (IC= 28.57-46.58 m/s,n = 8); discharge = 9.16 m’/s(IC = 3.32-21.65 m 5 /s). Aver¬ 
ages of physicochemical variables were: temperature = 28.03°C (1C = 27.02- 29.04°C, n = 30), pH = 
8.47 (IC = 8.33-8.61, n = 30), conductivity = 666.83 pS/cm (IC = 640.70 - 692.96 pS/cm, n = 30) and 
oxygen = 7.78 ppm (IC = 7.05 - 8.51 ppm, n = 30). 


t03*3ff 103*00' 102*30* 



MATERIAL AND METHODS 

COLLECTING. - Larvae were collected twice in each season (8 trips in total) at shores, and in 
riffles (at mid-channel) and eddies. Usually, sampling was done at the end of the first third and at the 
beginning of the last third of each season. We used an aquatic D-frame net, and samples were pre¬ 
served in 96% alcohol with one replacement before 24 h. A stereomicroscope was used to separate 
and quantify all larvae to the species level. In addition to the larval sampling, physicochemical vari¬ 
ables such as pH, dissolved oxygen, temperature and conductivity were recorded for each sample. 
Depth, width, and current velocity were measured and discharge was then calculated. The gradient 
(slope) was measured at seven points of the 500 m long sampling transect according to RESH et al. 
(1996). 

DIVERSITY MEASURES AND DATA ANALYSIS. — Richness and composition, as well as 
classical diversity measures such as the Shannon’s diversity index (IT), the Simpson’s index ( D ), the 
Margalef’s richness index ( R ), and the Pielou’s evenness index ( J ) were used in order to describe the 
Odonata assemblages by season, strata and as a whole (MORENO, 2001). Also, the Renyi’s diversity 
profiles (TOTHMERESZ, 1995,1998; JAKAB, 2002) were used for comparing diversity, as proposed 
by SOUTHWOOD & HENDERSON (2000). In this method when the value of the scale used as a 
parameter is low, the method is extremely sensitive to the presence of rare species. As the value of the 
scale increases, diversity is less sensitive to rare species. At a high value, the method is sensitive only 
to common species. The result of this scale-dependent characterization of diversity can be used in 
a graphical form to visualize the diversity relations of assemblages. This curve is usually called ‘the 




290 


J.A. Gomez-Anaya & R. Novelo-Gutierrez 


diversity profile of the assemblage’. It is important to stress that curves of two diversity profiles may 
intersect. For two communities, the intersection of the diversity profiles means that one of the com¬ 
munities is more diverse for rare species, while the other one is more diverse for common species. The 
Species Diversity & Richness package v. 3.0 was employed to generate the Renyi’s values, exporting 
them to an Excel spreadsheet to show them graphically. 

THEORETICAL RICHNESS. — An estimate of the theoretical richness using non paramet¬ 
ric estimators Chao2, Bootstrap, and upper limit of Mao Tau was carried out, using Estimates 8.0 
(COLWELL, 2006). Additionally, parametric methods as richness estimators by extrapolation were 
also used, which apply the observed curve of species accumulation for modeling the addition of new 
species in relation to the sampling effort (PALMER, 1990; SOBERON & LLORENTE, 1993). The 
Clench’s (CLENCH, 1979) and Linear dependence models were applied, as exemplified by JIMEN¬ 
EZ-VALVERDE & HORTAL (2000). Likewise, the slope on the cumulative species curve was used 
to assess the completeness of assemblages (HORTAL & LOBO, 2005). The slopes were obtained by 
means of the first derivative of the Clench’s and Linear dependence functions (NOVELO-GUTIER¬ 
REZ & GOMEZ-ANAYA, 2009). 

CLUSTER ANALYSIS. — Beta diversity was assessed by methods of classification. A Cluster 
Analysis (CA) on a Bray-Curtis (BC) similarity matrix [(1-W) where W = BC dissimilarity)] and the 
Unweighted Pair Group Method with Arithmetic mean (UPGMA) were used to explore the faunal 
relationships among seasons and strata. This analysis was performed using PC-ORD ver 4.5 (Mc- 
CUNE & GRACE, 2002). 

CANONICAL CORRESPONDENCE ANALYSIS. - The Canonical Correspondence Anal¬ 
ysis (CCA), a direct ordination method, was used to relate species abundance with environmental 
variables (TER BRAAK & SMILAUER, 1998). The number of environmental variables was then 
reduced using the automatic forward selection option in the CANOCO 4.5 program. The statistical 
significance of the relationship between the species and the set of environmental variables was tested 
by a Monte Carlo permutation test, using an F-ratio of the sum of all eigenvalues as the statistical 
test (TER BRAAK and PRENTICE, 1988). 

RESULTS 

SEASONAL PHYSICOCHEMICAL VARIATION 

Temperature, pH and conductivity had higher averages in summer, while ox¬ 
ygen was highest in spring. In fact, oxygen decreases gradually from spring to 
winter (Tab. I). All pH values were slightly alkaline. Oxygen levels were very low 
in winter when both abundance and species richness were highest. 

LARVAL RICHNESS AND COMPOSITION 

A total of 3,278 Odonata larvae belonging to 28 species (12 Zygoptera and 16 
Anisoptera), 16 genera and six families were collected (Tab. II). Erpetogomphus 
elaps (50.21%) was the dominant species; other numerically important species 
were Brechmorhoga praecox (14.16%) and Phyllogomphoides luisi (6.72%) (Fig. 
2). A further 61.54% of all species occurred in low abundance (<1%) and were 
considered rare. 



Table I 

Averages and 95% confidence intervals for the physichemical variables of Rio Pinolapa 
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SEASONAL NUMERICAL DOMINANCE 

During spring Erpetogomphus elaps (38.73%), Argia oe- 
nea (17.40%), and Progomphus marcelae (14.95%) domi¬ 
nated numerically. Only seven species were detected in 
summer (four of which were gomphids), with very low 
total and relative abundances. No coenagrionid was re¬ 
corded in this season. E. elaps (44.15%) and B. praecox 
(28.39%) dominated in autumn. Finally, during winter 
E. elaps (&).W/o) dominated the assemblage. Argia tezpi 
(8.74%) and Phyllogomphoides. luisi (6.52%) were present 
as codominant species together with E. elaps. It is in¬ 
teresting to note that while E. elaps clearly dominated 
throughout the year in the larval stage, the imagoes were 
not very commonly encountered. 

SEASONAL ASSEMBLAGES 

Table III and Figure 3 show the seasonal ecological pa¬ 
rameters. The smallest number of species was recorded 
in summer, spring and autumn were intermediate, and 
the highest number of species was found in winter. The 
abundance pattern seems to follow the richness one, be¬ 
ing higher in winter and lower in summer. Shannon’s 
diversity index (//’) seems to vary little throughout the 
four seasons, although it was a little higher in spring. 
Dominance ( D ) was higher in winter, mainly due to the 
great abundance of E elaps. Renyi's diversity profiles are 
shown in Figure 4. The summer pattern was a straight 
line. It showed the minimum number of species for a = 
0 (the basic structure of assemblages), but for values up 
to 2 (Simpson index) the pattern showed that summer 
diversity was the highest. 

SIMILARITY 

The Odonata larval assemblages from autumn and 
winter were the most similar, mainly due to their high 
and similar abundance (Fig. 5), nevertheless, they were 
quite different in richness sharing a high number of spe¬ 
cies (13). Some species like Argia oculata, Erpetogom- 
phus cophias and Paltothemis lineatipes were exclusively 



292 


J.A. Gomez-Anaya & R. Novelo-Gutierrez 


Table II 

Richness and composition of seasonal Odonata larval assemblages from Rio Pinolapa 


Taxa 

Key 

Spring 

% 

Summer % 

Autumn 

% 

Winter 

% 

Total 

% 

Number of individuals 


410 

12.45 20 

0.61 

1452 

44.32 

1396 

42.61 

3278 

Number of species 


15 

55.56 7 

25.93 

16 

59.26 

22 

81.48 

100 

ZYGOPTERA 

Calopterygidae 

Hetaerina 

americana 

Hearn 



4 20.0 

22 

1.52 

7 

0.50 

33 

1.01 

Platystictidae 

Palaemnema 

domina 

Pado 

23 

5.64 


16 

1.10 

54 

3.87 

93 

2.84 

Coenagrionidae 
Argia funcki 

Arfu 

7 

1.72 

. . 

. 

. 

. 

. 

7 

0.21 

A. oculata 

A roc 

- 

- 

- 

3 

0.21 

- 

- 

3 

0.09 

A. oenea 

A roe 

71 

17.40 

- 

14 

0.96 

53 

3.80 

138 

4.21 

A. pollens 

Arpa 

1 

0.25 

- 

- 

- 

8 

0.57 

9 

0.27 

A. pul la 

Arpu 

- 

- 

- 

- 

- 

6 

0.43 

6 

0.18 

A. tezpi 

Arte 

1 

0.25 

- 

20 

1.38 

122 

8.74 

143 

4.37 

Enallagma 

novaehispaniae 

Enno 

5 

1.23 

. . 

. 

. 

2 

0.14 

7 

0.21 

E. semicirculare 

Ense 

- 

- 

- 

- 

- 

3 

0.21 

3 

0.09 

Telebasis salva 

Tesa 

2 

0.49 

- 

- 

- 

- 

- 

2 

0.06 

Protoneuridae 

Protoneura cara 

Prca 






1 

0.07 

1 

0.03 

ANISOPTERA 

Gomphidae 

Erpetogomphus 

bothrops 

Erbo 






1 

0.07 

1 

0.03 

E. cophias 

Erco 

- 

- 

- 

1 

0.07 

- 

- 

1 

0.03 

E. elaps 

Erel 

158 

38.73 

7 35.0 

641 

44.15 

839 

60.10 

1645 

50.21 

Progomphus 

clendoni 

Prcl 

13 

3.19 

2 10.0 

9 

0.62 

8 

0.57 

32 

0.98 

P lambertoi 

Prla 

1 

0.25 

- 

91 

6.27 

1 

0.07 

93 

2.84 

P. marcelae 

Prma 

61 

14.95 

- 

8 

0.55 

65 

4.66 

134 

4.09 

Phyllogomphoides 
luisi Phlu 

21 

5.15 

2 10.0 

106 

7.30 

91 

6.52 

220 

6.72 

P. pacificus 

Phpa 

- 

- 

2 10.0 

9 

0.62 

8 

0.57 

19 

0.58 

Libellulidae 

Dythemis 

nigrescens 

Dyni 






7 

0.50 

7 

0.21 

Brechmorhoga 

praecox 

Brpr 

5 

1.23 

2 10.0 

412 

28.37 

45 

3.22 

464 

14.16 

Erythrodiplax sp. Ersp 

- 

0.00 

- 

- 

- 

7 

0.50 

7 

0.21 

Macrothemis 

inacuta 

Main 

7 

1.72 






7 

0.21 

M. pseudimitans 

Maps 

28 

6.86 

1 5.0 

83 

5.72 

60 

4.30 

172 

5.25 

Paltothemis 

lineatipes 

Pali 



. 

3 

0.21 

. 

. 

3 

0.09 

Perithemis 

domitia 

Pedm 






2 

0.14 

2 

0.06 

Pseudoleon 

superbus 

Pssu 

6 

1.47 

- 

14 

0.96 

6 

0.43 

26 

0.79 
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Fig. 2. Relative abundance of Odonata species from Rio Pinolapa. Key to species in Table II. 


recorded in autumn, while others as Argia pulla, Enallagma semicirculare, Fro- 
toneura cara, Dythemis nigrescens and Erythrodiplax sp., were only found in win¬ 
ter. Summer contains an assemblage of very few species. 

DIVERSITY AND ABUNDANCE BY STRATA 

Table IV shows the ecological parameters of the Odonata larval assemblages 
by strata of RP. Number of species in the shores was almost twice as much of 
that of the riffles (middle-channel), and more than twice as much of that of ed¬ 
dies. However, the abundance was higher in the middle-channel of the water body. 
The diversity H’ was higher in the shores while the dominance did in eddies. 

The major amount of larvae in all strata was Erpetogomphus elaps, however, 
in the riffles we found a higher proportion of Brechmorhoga praecox (24.42%). 

When diversity is compared and ordered using the Renyi’s profiles (Fig. 6), it 


Table III 

Seasonality of the ecological parameters of Odonata larval assemblages at Rio Pinolapa 



Spring 

Summer 

Autumn 

Winter 

Total 

Number of species 

15 

7 

16 

22 

28 

Number of specimens 

401 

20 

1452 

1396 

3269 

Simpson (Z>) 

0.22 

0.21 

0.29 

0.38 

0.29 

Shannon-Weaver (FT) 

1.89 

1.76 

1.63 

1.61 

1.85 

Margalef richness ( R ) 

2.34 

2.00 

2.06 

2.90 

3.09 

Pielou evenness ( J) 

0.70 

0.90 

0.59 

0.52 

0.57 
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Fig. 3. Ecological parameters per season of the Odonata larval assemblage at Rio Pinolapa. Collec¬ 
tions were made between March 2005 and January 2006. 


follows the gradient shores>riffles>eddies. This fact confirms that shores contain 
the highest diversity of Odonata larvae in RP. 

THEORETICAL RICHNESS ESTIMATION 


non-parametric MODELS. — The cumulative species curves generated by 
non parametric estimators Chao2, and Bootstrap are shown in Figure 7. The es¬ 
timated number of species was 41.4, and 32.4 species, respectively, which gave a 
sampling efficiency of 67.6%, and 86.3, respectively. These richness estimators 
indicate a lack of register from 4 to 13 species. The estimated number of species 
using the Mao Tau upper limit of class interval was 34.7, which means that still 
should be added to the list 6-7 species, being the efficiency of the total sampling 
effort of 80%. The number of species with a single individual (singletons) was 3, 
with two individuals (doubletons) was 2, the number of unique species was 11, 
and for duplicated ones was 3. 


PARAMETRIC MODELS. - 
Figure 8 shows cumulative 
species curves generated by 
the Clench’s and Linear de- 


Table IV 

Ecological parameters of the Odonata larval assemblages 
by stratum at Rio Pinolapa 


pendence functions. The first 
function predicted 30.67 spe¬ 
cies and explained 98% of 
data variation, while the sec¬ 
ond one predicted 25.74 spe¬ 
cies and explained 95% of 
data variation. The Clench’s 
model indicates that there are 
3 species to be registered yet, 
and the linear dependence in- 


Index/stratum 

Shores 

Riffles 

Eddies 

Number of species 

25 

14 

10 

Number of specimens* 

1062 

1713 

323 

Simpson ( D ) 

0.27 

0.32 

0.47 

Shannon-Weaver (If) 

1.93 

1.54 

1.18 

Margalef’s richness (R) 

3.44 

1.61 

1.56 

Pielou’s (J) eveness 

0.60 

0.60 

0.51 


* Some samples were omitted because of insufficient field 
data. 
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Fig. 4. Renyi’s diversity profiles for the four seasonal Odonata larval assemblages from Rio Pinola¬ 
pa. 


dicates the list has been completed. Based on the explained variance ( R 2 , deter¬ 
mination coefficient), the Clench’s estimation explained better the data variation; 
its prediction is considered further. Finally, slopes for both curves were 0.10 for 
Clench’s function and 0.04 for Linear dependence. 

THE SPECIES-ENVIRONMENT RELATIONSHIPS 

The results of the CCA were globally significant (trace = 0.937, F = 1.59, p 
<0.05, Tab. V). The first three axes offered a good solution to the ordination of 
the physicochemical variables and abundance of species, since from the total vari¬ 
ability in the data (inertia = 3.701), it was possible to explain 87.9% by means of 
these group of axes. The significance test of the first canonical axis showed it was 
significant (eigenvalue = 0.292, F = 2.822, p <0.05). 


Distance (Objective Function) 

6 IE-02 2.3E-01 4E-01 5.7E-01 7.4E-01 

Information Remaining (%) 

100 75 50 25 0 

Spring - 

Autumn i_ 

Winter ■ 

Summer- 

Fig. 5. Dendrogram showing the faunal relationships among the year seasons assemblages of odo- 
nate larvae. Based on a Bray-Curtis similarity matrix and the unweighted pair-group arithmetic av¬ 
eraging rUPGMA). 
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Table V 

Results of the canonical correspondence analysis (CCA) of log-transformed Odonata larvae 
abundance as a function of their environmental variables 


Axes 

1 

2 

3 

4 

Total inertia 

Eigenvalues 

0.29 

0.25 

0.18 

0.09 

3.705 

Species-environment correlations 
Cumulative percentage variance 

0.82 

0.85 

0.76 

0.66 


of species data 

7.9 

14.7 

19.5 

22 


of species-environment relation 

Sum of all eigenvalues 

Sum of all canonical eigenvalues 

31.2 

58 

76.9 

86.8 

3.705 

0.937 


The first axis was the most important, explaining 31.2% of variance, and it was 
also the most strongly correlated with pH, riffles and autumn. The second axis 
explained 26.8% of variance, and it correlated strongly with conductivity. The 
third and fourth axes explained only 16.7% and 12.1% of variance, respectively, 
and were not considered further. 

When the distribution of the species in the eight collections and three strata 
is analyzed together with the CCA of the Figure 9, it is posible to make the fol¬ 
lowing precisions: Paltothemis lineatipes, Erpetogomphus cophias, Progomphus 


3.5 


3 



0 

0.25 0.5 0.75 1 1.25 1.5 1.75 2 2.25 2.5 2.75 3 3.25 3.5 3.75 4 
Scale parameter 

Fig. 6. Renyi’s diversity profiles for three strata assemblages of Odonata larvae from Rio Pinolapa. 
Profiles differ mainly at their basic level of structure, the number of species. Profiles never cross in 
a range from 1 to 4. Values of Renyi when a = 4 were: shores profile = 0.982, riffles profile = 0.969 
and eddies profile = 0.535. 
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lambertoi, and Hetaerina Americana do associate well with the autumn. Particu¬ 
larly, P. lambertoi was registered almost in 100% (91 out of 92 specimens) from 
the shores at the beginning of the autumn. Argia oenea and Palaemnema domina 
were almost invariably registered from the riffles and at late winter. Only seven 
larvae of Macrothemis inequiunguis were registered from late spring on shores. 

DISCUSSION 

GENERAL ASPECTS 

A total of 28 Odonata species (12 Zygoptera and 16 Anisoptera) were found 
and 3,276 larvae were identified in this work. The size of the Odonata larval as¬ 
semblage from RP is similar to other Mexican water bodies reported (ALONSO- 
EGUIALIS, 2004; GOMEZ-ANAYA et al., 2000; NOVELO-GUTIERREZ & 
GONZALEZ-SORIANO, 1991; NOVELO-GUTIERREZ et al., 2002; BOND 
et al., 2006). Usually, the most speciose families are Libellulidae and Coenagrio- 
nidae. In this case Coenagrionidae, Gomphidae and Libellulidae were best rep¬ 
resented. The single genus Argia contributed 50% of Coenagrionidae; in Gom¬ 
phidae three genera contributed more or less equally, and in Libellulidae seven 
genera contributed to the diversity. 

DIVERSITY AND ABUNDANCE 

The highest and lowest species richness was found in winter and summer, re¬ 
spectively. The greatest abundance of individuals was recorded in autumn and 

45 n 



1 3 5 7 9 11 13 15 17 19 21 23 25 27 29 31 33 35 37 39 41 

Sampling effort 


Fig. 7. Cumulative species curves generated by non-parametric estimators Mao Tau (S ob5 ), Chao2, 
and Bootstrap for the Odonata larval assemblage from Rio Pinolapa. 
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1 2 3 4 5 6 7 8 9101112131415161718192021222324252627282930313233343536373839404142 
Sampling effort 

Fig. 8. Cumulative species curves generated by the Clench’s function 
(a = 3.68, b = 0.12, asymptote = 30.67, R = 0.98), and linear depend¬ 
ence function (a = 2.63, b = 0.102, asymptote = 25.74, R = 0.95), for 
the Odonata larval assemblage from Rio Pinolapa. The slope was esti¬ 
mated by the first derivative of Clench’s function [a/(l+b*n) 2 ] and the 
first derivative of Linear dependence function (a*exp(-b*n)]. 


winter. This could be 
due to the emergence 
pattern of the sea¬ 
sonal species. Most 
of them emerge dur¬ 
ing the dry season (at 
the end of the win¬ 
ter and through the 
spring). Thus, when 
the rains come in the 
summer, very few spe¬ 
cies are present as lar¬ 
vae. Summer is the 
season with a great 
reproductive activity, 
so that when autumn 


comes, there are large populations of larvae of the majority of odonate species. 
Later, as time goes on, abundance decreases through mortality until the lower 
abundances in spring and then, with emergence of the adults, even lower in sum¬ 
mer. Among strata, the highest species richness was found at shoreline and the 
lowest in eddies; while the highest abundance occurred in riffles and the lowest 
in eddies. Eddies were the most simplified stratum, usually with a muddy bottom 
and some decaying leaves, and lacking any kind of aquatic plant; this could result 
in the lowest richness and abundance. Richness distribution pattern was differ¬ 
ent to that of abundance in RR Thus, 55.39% of all abundance was recorded in 
riffles, 34.20% in shorelines and 10.40% in eddies. The most abundant species in 
all three strata was Erpetogomphus elaps, with 53.79%, 32.44%, and 13.77%, re¬ 
spectively. This species is, apparently, the best adapted to different conditions in 
time and space in RP, despite having a restricted period as adults (late summer 
to early autumn). 


RENYI'S DIVERSITY PROFILES 


Diversity ordering can be performed by using a diversity index family (TOTH- 
MERESZ, 1995), as proposed by SOUTHWOOD & HENDERSON (2000). The 
Renyi’s diversity profiles method has been poorly used in comparing and order¬ 
ing diversity of aquatic macroinvertebrate assemblages (SIPKAY et al. 2007) and 
even less used with Odonata assemblages (JAKAB et al., 2002). In this method a 
scale parameter (alpha) is related to the abundance-dominance structure of the 
community. At different values of the scale the function is sensitive to rare, com¬ 
mon, and intermediate-abundance species. When comparing two Renyi diversity 
profiles, if they do not cross each other, it means that the upper profile represents 
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a more diverse assemblage under any common measure of diversity; but if these 
profiles do cross each other once, it means that one assemblage is more diverse 
when rare species are weighted more heavily (low alpha), and the other assem¬ 
blage is more diverse when common species are weighted more heavily (high al¬ 
pha). When the diversity profiles cross each other, the communities cannot be 
ordered according to their diversity, because one of them is more diverse for rare 
species, the other for common ones. Renyi’s diversity is sensitive to rare species 
for small values of the scale, whereas it is sensitive to abundant species for larger 
values of the scale. Seasonal diversity profiles crossed once, mainly for a<l val¬ 
ues (Fig. 4). It means that these assemblages are mainly different in rare species. 
Strata profiles do not cross; however, shore and riffles profiles were similar in 
abundant and frequent species (a>2), and, definitively, eddies profile was the least 
diverse. 

Species distribution must meet the ecological requirements of all stages in the 
life cycle (CORBET, 1999). Although the imagoes have more mobility than imma¬ 
ture, adequate substrates must be present for the larval emergence (rocks, vegeta¬ 
tion, twigs, etc.). In consequence, the disturbance of the original conditions (e.g. 
riparian vegetation) will affect richness (SMITH et al., 2007). In RP the major 
diversity in shores can be explained because of major microhabitat heterogene¬ 
ity. In shores we observed different size of rocks, sand, algae, plants, roots, mud, 
leaf packs, detritus, and different combinations of these substrates. Number of 
substrates in riffles and eddies or combinations of them were limited. In shores 
we found 10 Zygoptera species (1 calopterygid, 1 platystictid, 1 protoneurid, 7 
coenagrionids), and 15 Anisoptera (7 gomphids and 8 libellulids). All Zygoptera 
species are endophitic in oviposition and they can find the needed substrates in 
shores. Gomphids and libellulids females lack ovipositor and they release their 
eggs on the water surface near the shores. 

Most of the riffle samples of RP come from shallow and smooth flowing wa¬ 
ter places. In riffles we found 14 species (5 Zygoptera and 9 Anisoptera). Brech- 
morhoga praecox was the most abundant species in this stratum, as have been re¬ 
ported to different species of this genus (CORBET, 1999), and particularly for 
this species (BOND et al., 2006). Most of B. praecox larvae were caught, main¬ 
ly, at the end of autumn. They exhibit different disruptive color patterns, being 
cryptic at sand and fine gravel bottom. Must of Coenagrionidae species were re¬ 
corded from shores, except Argia oenea and A. tezpi which were more abundant 
in riffles. 


CUMULATIVE CURVES OF SPECIES 


The non-parametric estimation achieved by Bootstrap agrees well with the par¬ 
ametric estimation made by Clench function, which predicted 31 species (3 spe¬ 
cies more to be added). For parametric estimations, we chose initially the model 
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that better fitted data using the coefficient of determination (R 2 ). This model 
was the Clench’s function. Then, we used the slope of the curve evaluated at the 
maximum sampling effort by the first derived of both functions, as an approach 
of completeness of the list, according to JIMENEZ-VALVERDE & HORTAL 
(2003), and HORTAL & LOBO (2005). Although the slope was smaller with the 
Linear dependence (0.04), we believe this function underestimates the number 
of species (since it predicted fewer species than the recorded ones). We think few 
species could exist yet to be added to the list for the studied section of the river. 
The slope of the first derived of Clench (0.10) is relatively small, supporting the 
idea that the list is practically complete. 

CCA 


Several authors have explored the relationship between some environmental 
factors and odonate faunas using multivariate techniques (SAMWAYS, 2003 
SCHINDLER et al., 2003; HOFFMANN & MASON, 2005; OPPEL, 2005 
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LEN, 2008; SATO & 
RIDDIFORD, 2008; 
HAMASAKI et al., 
2009). In this work the 
CCA showed signifi¬ 
cant correlation be¬ 
tween both environ¬ 
mental and species 
0.8 matrices. Additional¬ 
ly, the CCA biplot al¬ 
lowed visualization of 
some particular facts 
and the establish¬ 
ment of some hypoth¬ 
eses on species-species 
and environment-spe¬ 
cies relationships. For 
example, Palaemnema 
domina and Argia oe- 
nea are well associat- 


Fig. 9. Biplot of CCA ordination showing environmental variables (ar- e£ j riffles Stratum 
rows) most strongly correlated with axes CC1 and CC2 and species in indicad S]milar ^ 
triangles. Seasons and strata were included in CCA as dummy variables. . ~ ^ 

In terms of predicting larval assemblage composition, important envi- logics' requirements, 
ronmental variables have longer arrows than less important ones. P domina is usually 
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found under mid-sized rocks (8-10 cm diam.) in riffles, as is A. oenea, although 
this last species is found also commonly among rough gravel. Erpetogomphus elaps 
is close to the origin on Figure 9, meaning it is the most ubiquitous species in RP. 
Apparently, the larvae are well adapted to both erosional and depositional envi¬ 
ronments. Moreover, it was found in every collection and through all seasons. 

FINAL CONSIDERATIONS 

Contrary to species lists from lentic water bodies (lakes, lagoons, ponds) which 
can be considered complete lists because one can sample the whole water body, 
the 28 recorded species in the RP survey represent just a part of a bigger assem¬ 
blage. This bigger assemblage extends up and down stream changing in richness, 
composition, and abundance with changes in river conditions. From this point 
of view, the assemblage here described represents a relatively local measure of 
this group of insects, and also a measure of the conservation status of the river. 
It is possible that for many species with low abundance and distribution in space 
(strata) and time (seasons), the best conditions for their reproduction could be 
up or down stream. In this bigger assemblage the floristic, climatic, altitudinal 
and microhabitat changes should be considered, because the more the changes 
in the river conditions, the larger the species list should be. 
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